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Abstract 

This paper studies particle propagation in a one-dimensional inho- 
mogeneous medium where the laws of motion are generated by chaotic 
and deterministic local maps. Assuming that the particle's initial loca- 
tion is random and uniformly distributed, this dynamical system can 
be reduced to a random walk in a one-dimensional inhomogeneous en- 
vironment with a forbidden direction. Our main result is a local limit 
theorem which explains in detail why, in the long run, the random 
walk's probability mass function does not converge to a Gaussian den- 
sity, although the corresponding limiting distribution over a coarser 
diffusive space scale is Gaussian. 

Keywords: random walk in random environment, quenched random walk, 
random media, local limit theorem, extended dynamical system 
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1 Introduction 

1.1 A chaotic dynamical system 

This paper studies a particle moving in a continuous inhomogeneous medium 
which is composed of a linear chain of cells modeled by the unit intervals 
[k, A: -|- 1) of the positive real line. Each interval [k, A: -|- 1) is assigned a label 
LOk and a map C/^j. which determines the dynamics of the particle as long as 
the particle remains in the interval. The sequence of labels u = {oJkjk&z+i 
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called an environment, is assumed to be either nonrandom, or a realization 
of a random sequence that is frozen during the particle's lifetime. 

We are interested in the case in which the local dynamical rules [7^^. are 
chaotic in the sense that the distance between two initially nearby particles 
grows at an exponential rate. More concretely, we shall focus on a model 
where a particle located at x„ G [k,k + 1) at time n jumps to Xn+i = 
k + Ui^^{xn — k). Here oj^ G (0, 1) and C/^j. is the piecewise affine map from 
[0, 1) onto [0, 2) such that U^^ [0, 1 - Wfc) = [0, 1) and C^j^ [1 - Uk, 1) = [1, 2). 
The dynamical system generated by the local rules is compactly expressed 
by Xn+i = U^ixn), where the global map on the positive real line is 
defined by 

U^{x) = [x] + U^^^^{x-[x]), (1.1) 
and [x\ denotes the integral part of x, see Figure 1. 




Figure 1: The piecewise affinely on the interval [k,k + 1). 



The above model belongs to the realm of extended dynamical systems, a 
somewhat vaguely defined yet highly active field of research (e.g. Chazottes 
and Fernandez [9]). Telltale characteristics of such systems are a noncom- 
pact or high-dimensional phase space and the lack of relevant finite invariant 
measures. Our principal motivation is to study the impact of environment 
inhomogeneities on the long-term behavior of extended dynamical systems. 
Concrete models include neural oscillator networks (Lin, Shea-Brown, and 
Young and the Lorentz gas with randomly placed scatterers (Chernov 
and Dolgopyat [lj]; Cristadoro, Lenci, and Seri [12]), to name a few. In this 
paper, we shall restrict the analysis to the affine dynamical model in (1.1), 
to keep the presentation simple and clear. 

1.2 Random initial data 

Because the local maps Uuj^. are chaotic, predicting the particle's future loca- 
tion with any useful accuracy over any reasonably long time horizon would 
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require a precise knowledge of its initial position — a sheer impossibility in 
practice. Therefore, it is natural to take the statistical point of view and 
study the stochastic process defined by 



To analyze the time evolution of the above process, we must impose 
some regularity conditions on the environment. In particular, those condi- 
tions guarantee ballistic motion, and one might guess that the distribution 
approaches Gaussian in the long run. To test this hypothesis, we have plot- 
ted in Figure 2 numerically computed histograms of Xn at time n = 2^^ in 
two frozen environments, using the intervals [k,k + 1) as bins. Rather sur- 
prisingly, the histograms do not appear Gaussian. A similar phenomenon 
was recently observed by Simula and Stenlund [33, 34]. 

1.3 Summary of main results 

The main contribution of the paper is to explain the emergence of the his- 
tograms in Figure 2. This is accomplished by first reducing the continuum 
dynamical system to a unidirectional random walk on the integers (Theo- 
rem 2.1), and then deriving a local limit law (Theorem 2.4) that completely 
explains the behavior observed in Figure 2. As a byproduct, we also obtain 
a law of large numbers (Theorem 2.3) and a central limit law (Theorem 2.5) 
for the walk. These limit laws are valid for all frozen environments — ran- 
dom or nonrandom — which satisfy certain statistical regularity properties. 
We also devote a separate section to the analysis of random environments, 
where we show (Theorem 2.7) that the three aforementioned limit laws are 
valid for almost all realizations of a stationary random environment under 
suitable moment and mixing conditions. Because the random walk is uni- 
directional (it never steps backwards), limit theorems for its hitting times 
are immediate consequences of classical limit laws for independent random 
variables. Translating the limit theorems of the hitting times into limit the- 
orems of the walk location form the main task in proving the results; see 
Section 3.3 for a general outline of the proofs. 

1.4 Related work 

We shall discuss here only literature most closely related to transient one- 
dimensional random walks in quenched random environments; for a broad 
picture of the theory of random walks in random environments, see e.g. 
Bolthausen and Sznitman [6] , Sznitman [38] , and Zeitouni [ ] . Laws of large 
numbers and averaged central limit theorems for random walks in random 
environments have been known already for a long time (e.g. Solomon [35]; 
Kesten, Kozlov, and Spitzer [_ j]), whereas the literature on quenched central 
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Figure 2: Histograms of x„, (black dots) at time n = 2^^ in two frozen 
environments uj. Top: lj is a realization of a Markov chain with values 
in {5,1} and P(a;fc+i = w^lcufc) = |. Bottom: uj is nonrandom, with 
= ^ + ^ sin/c. The green circles are obtained by modulating the blue 
Gaussian density by the factor uj'^^ / fi appearing in Theorem 2.4. 
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limit theorems is more recent. Buffet and Hannigan [S] proved a quenched 
central limit theorem for a pure birth process in an independently scattered 
random environment under moment conditions later relaxed by Horvath 
and Shao [21, 22]; this model is a direct continuous-time analogue of the 
random walk studied here. Quenched central limit theorems for more gen- 
eral one-dimensional transient random walks were proved only very recently, 
independently by Goldsheid [ ] and Peterson [30] (see also Alili [1] for a 
result concerning a special quasiperiodic environment). Rassoul-Agha and 
Seppalainen ["52] obtained a similar result for multidimensional random walks 
with a forbidden direction, which in the one-dimensional case corresponds to 
the unidirectional walk analyzed in this paper. Dolgopyat, Keller, and Liv- 
erani [14] have obtained a quenched central limit theorem for environments 
changing in time and space. 

Our approach differs from most earlier works in that we separately an- 
alyze the two degrees of randomness involved in random walks in quenched 
random environments. In the first part, we extract a set of statistical regu- 
larity properties for a given environment that are sufficient for proving the 
limit laws, while treating the environment as nonrandom. In the second 
part, we derive conditions for the probability distribution of the random en- 
vironment that yield almost surely regular realizations. A key result for the 
second part is a law of large numbers for the moving averages of a stationary 
sequence (Lemma 4.2), which is proved with the help of Peligrad's exten- 
sion [29] of the Baum-Katz theorem [3] (see Bingham [ I] for a nice survey). 
The major advantage of this approach is a clarified picture on how differ- 
ent sources of randomness affect the random walk's behavior in quenched 
random environments. 

Local limit theorems for random walks in homogeneous environments 
(e.g. Spitzer [ ]) can be proved as simple consequences of Gnedenko's clas- 
sical theorem (e.g. [16, Chapter 9] or [1-5, Section 3.5]), because such walks 
are just sums of independent random variables. A generalization for a pe- 
riodic environment was derived by Takenami [■' ], and for a periodic graph 
recently by Kazami and Uchiyama [_ .]. In contrast, local limit theorems 
for random walks in aperiodic inhomogeneous environments appear nonex- 
istent. To the best of our knowledge. Theorem 2.4 and Theorem 2.7 are the 
first local limit laws concerning transient random walks in aperiodic nonran- 
dom or quenched random environments. Although our analysis is restricted 
to a very special instance of a random walk, the model is still rich enough 
to capture several interesting phenomena, such as the need for nonlinear 
centering and a non-Gaussian modulating factor, and we believe that the 
results could serve as useful benchmarks when testing hypotheses concerning 
more general random walks. 

Regarding extended dynamical systems, we have found two earlier local 
limit theorems, both corresponding to homogeneous environments. Szasz 
and Varjii [37] have considered Lorentz processes with periodic configura- 
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tions of scatterers, while Bardet, Gouezel, and Keller [2] study rather dif- 
ferent type of systems: small (possibly inhomogeneous) perturbations of 
weakly coupled, translation invariant, coupled map lattices; see Nagaev [28] 
and Guivarc'h [20] for some of the original techniques. Let us finally stress 
that for more classical, probability measure preserving, dynamical systems, 
various types of limit theorems have been proved for many decades. Yet 
such systems, too, continue to be studied vigorously, with important recent 
developments (e.g. Chazottes and Gouezel [!()]; Gouezel [18, 19]). 

1.5 Organization of the paper 

The rest of the paper is organized as follows. Section 2 presents the main 
results. The proofs for nonrandom environments are given in Section 3, 
and the proofs for quenched random environments in Section 4. Section 5 
concludes the paper, and Appendix A contains basic facts on generalized 
inverses of increasing sequences. 

2 Main results 

2.1 Representation as a random walk 

The distribution of the dynamical system (1.2) at any time instant can 
be completely characterized in terms of the following simple unidirectional 
random walk (discrete-time pure birth process) on the integers. Given an 
environment u; G (0, 1)^+ , let {Xn)neZ+ be a random walk in having the 
initial state Xq = and transitions 



We denote by P^j the distribution of the walk in the path space Z_,_^ . Note 
that if the environment a; is a realization of a random sequence, the process 
(Xn) can be identified as a random walk in a random environment, and the 
distribution P^^ is usually called the quenched law of the random walk. The 
expectation and variance with respect to P^j are denoted by and Var^^, 
respectively. When presenting general facts in probability theory, we write 
P and E for the measure and expectation. 

Theorem 2.1. For any environment lo, the value of the dynamical sys- 
tem (1.2) at any time instant n has the same distribution as Xn + R, where 
{Xn) is the random walk in Z+ defined by (2.1), and R is a uniformly dis- 
tributed random variable in [0, 1) independent of X^. 

Proof. The proof follows by induction and Lemma 2.2 below. □ 




A;, with probability 1 — ujk 
k + 1, with probability uj^. 



(2.1) 
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Lemma 2.2. Let X be a positive random variable such that (i) [X] and 
{X} = X — [X] are independent, and (ii) {X} is uniformly distributed in 
[0, 1). Then the same is true forU^X, and moreover, 



Pi[U^X] = l\[X] = k) 



1-Uk, ifl = k, 
UJk, ifl = k + l, 



(2.2) 



whenever P{[X] = k) > 0. 



Proof. Denote X = K + R, where X is a positive random integer indepen- 
dent of R, and R is uniformly distributed in [0, 1). Assume first that K = k 
for some nonrandom integer k. Then hli,jX = k + Ui^^R. Moreover, a simple 
calculation based on the definition of U^jj, shows that for all r G [0, 1), 



P{[U^X] = I, {U^X] < r) 



'{l-ujk)r, ifl = k, 

Ukv, if / = + 1, (2.3) 
0, else. 



Hence Ui^X satisfies (i), (ii), and (2.2) in the case where [X] is nonrandom. 
The general case follows by conditioning on X. □ 

2.2 Limit theorems for regular frozen environments 

In this section we shall analyze the random walk (X„) in a fixed, sufficiently 
regular environment, which may either be nonrandom, or a realization of 
a random sequence. More precisely, we shall in general assume that the 
environment ui G (0, 1)^+ is such that 



^k' = 0{k') 



for some < A < 1/2; and 



fc-i 



j=0 



k-l 



k~' ^(1 - ojj)ojf = a^ + o{k-\\og fc)-i/2), 

j=0 

for some constants /x > 1 and cr^ > 0. Moreover, we assume that 

k-l 

j=0 



(2.4) 

(2.5) 
(2.6) 



(2.7) 



and that the environmental moving averages satisfy 

o(A;V2-A) 



max 

j:\j\<ub{k) 



k+j-l 
£=k 



(2i 
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for all u > 0, where b{k) = (klogk)^^"^ . (In the special case with A = it 
suffices to use b{k) = /c^/^). Concrete examples of environments that satisfy 
the above regularity conditions shall be given in Section 2.3. 

The quantity cj^^ represents the mean sojourn time of the particle in the 
interval [k, k+1) (see Section 3.1 for more details). Therefore, the constant /i 
appearing in (2.5) may be interpreted as the inverse of the particle's traveling 
speed. The following result confirms this intuition. 

Theorem 2.3 (Law of large numbers). For any environment lo satisfy- 
ing (2.4) - (2.5), 

PUn-^Xn^ fi-^) = l. (2.9) 

The main result of the paper is the following limit theorem for the ran- 
dom walk {Xn) defined by (2.1), or alternatively (by virtue of Theorem 2.1), 
for the dynamical system defined by (1.2). 

Theorem 2.4 (Local limit theorem). For any environment ui satisfying (2.4 ) 
- (2.8), 

PUXn = k) = ^^^e-^^5i^ + o(n-V2), (2.10) 

uniformly with respect to k > 0, where / fi^ , and the centering factors 

k'^ are given by 

/fc^ = min ^ > : ^ > n V . (2.11) 

Two features in Theorem 2.4, which distinguish it from classical limit 
theorems, call for special attention. First, the centering factors fcj^ depend 
on the environment, and are in general nonlinear functions of n. Second, 
the modulating factor uj^^ / ^ in (2.10) causes the asymptotic shape of the 
probability mass function of X„ to be non-Gaussian. This modulating factor 
explains the behavior observed in Figure 2. 

In contrast, when looking at the probability distribution of the walk 
over a coarser diffusive space scale, the non-Gaussian modulating factor in 
Theorem 2.4 averages out asymptotically, and we end up with a standard 
Gaussian limiting distribution. 

Theorem 2.5 (Central limit theorem). For any environment uj satisfy- 
ing (2.4) - (2.8) for some constants fj, > 1 and o"^ > 0, 

P. ( <x]^^ r e-*V2 dt (2.12) 
V o-^n J V2vr J-00 

for all X E M, where o"^ = a'^/fj,^, and k'^ are given by (2.11). 
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Remark 2.6. The centering factor /cj^ may be identified as a generalized 
inverse (Appendix A) of the function k i— )• E^^T^-, where denotes the 
hitting time of the walk into site k (Section 3.1). A quick inspection of the 
proofs in Section 3 shows that Theorem 2.4 remains true if k^ is replaced 
by = k^ + o{'n}^'^^^), where A is as in (2.4); use the inequality |e~^ — 
e^y I — y\ for k <n and the proof of Lemma 3.6 for k > n. Moreover, 
Theorem 2.5 remains true if fej:^ is replaced by = k':^ + o(n-^/^). 

2.3 Limit theorems for quenched random environments 

In this section we assume that the environment a; is a realization of a sta- 
tionary random sequence in (0, 1)^+ , and denote by P its distribution on 
(0,1)^+. The expectation with respect to P is denoted by E. We shall 
assume that 

¥.{ujQ^f < oo for some g > 5. (2.13) 

To guarantee that the environmental averages converge to their mean values 
rapidly enough, we assume that 

j;0i/2(A;)<oo, (2.14) 
fe>i 

where the mixing coefficients (j){k) are defined by 

(/>(A;)=sup sup |P(5|A) -P(B)|, (2.15) 

and where = cr{uj,j < m) and = cr{u}j,j > m) (e.g. Bradley [7]). 

The following result summarizes three limit theorems for the quenched 
random walk in a stationary strongly mixing random environment. 

Theorem 2.7. The law of large numbers (2.9), the local limit law (2.10), 
and the central limit law (2.12) are valid with /i = E(^(^^ and = E(l — 
ujq)uJq for almost every realization of a stationary random environment 
satisfying (2.13) and (2.14)- 

Especially, the limit laws summarized by Theorem 2.7 hold in the fol- 
lowing cases: 

• Independently scattered stationary environments (environments where 
the site labels oJk are independent and identically distributed). 

• Uniformly ergodic environments as discussed in Goldsheid [17]. 

• Environments which are realizations of finite-state irreducible aperi- 
odic stationary Markov chains, or more general Markov chains satis- 
fying Doeblin's condition (e.g. [7]). 
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Although in many apphcations it is natural to assume that the environment 
is stationary, the limits of Theorem 2.7 remain valid under looser conditions, 
as is clear from the results of Section 2.2. 

Alternative versions of the law of large numbers and the central limit 
law in Theorem 2.7, where the centering factors fcj^ are replaced by E^^Xn, 
can be proved as consequences of Theorems 3.1 and 5.4 in Rassoul-Agha 
and Seppalainen [.32], if we additionally assume that 



This so-called nonnestling assumption is close in spirit to the uniform el- 
lipticity of nearest-neighbor random walks in random environments; in the 
context of our model it corresponds to the special case A = in (2.4). Al- 
though we believe that the local and central limit laws in Theorems 2.4, 2.5, 
and 2.7 remain generally valid also for the alternative centering n i— ?■ E^X^, 
we prefer to use the centering k'^ defined in (2.11), because these factors 
are easily computed from the environment. Analogous central limit laws for 
nearest-neighbor walks were recently independently found by Goldsheid [17] 
and Peterson [.')()]. 

If we were only interested in the law of large numbers (2.9), we could do 
with less assumptions in Theorem 2.7. For example, as our proof in Section 4 
shows, the moment condition (2.13) would only be needed for q > 2. The 
mixing assumption (2.14) could be relaxed as well, see for example Bingham 
[4] for more details. 

3 Proofs for regular nonrandom environments 

This section is devoted to proving Theorems 2.3-2.5 for the random walk 
(Xn) in an environment that satisfies the regularity assumptions (2.4) - 
(2.8). The environment u shall be fixed once and for all during the whole 
section — here we do not care whether it is a realization of a random se- 
quence or not. 

The section is organized as follows. Section 3.1 describes some prelim- 
inaries on the hitting times of the walk, and Section 3.2 gives the proof of 
the law of large numbers. The proof of the local limit theorem is split into 
Sections 3.3-3.7, and the proof of the central limit theorem is in Section 3.8. 

3.1 Hitting times of the walk 

We denote the hitting time of (X„) into site /c by = min{n > : Xn = k}, 
and the sojourn time at site k hy = T^+i — T^. Because the walk never 
moves backwards, the equivalence 




for some 5 > 0. 



Xn = k if and only if <n < T^+i 



(3.1) 
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is valid for all k and n. Moreover, the sojourn times are independent, and 
has a geometric distribution on {1, 2, ... } with success probability w^. Hence 
the mean and the variance of Tk are given by -E'^j'^fc = and Var^(Tfc) = 
(1 — Wfc)^^^, respectively. The mean and the variance of are denoted by 
Hk = E^Tk and cr^ = \ai^{Tk), so that 

Hk = X]^7^ '^k = ~ ^i)^^7^- (3.2) 

i=o i=o 

The following result transforms the realization-by-realization relationship 
(3.1) into one concerning the probability mass functions. 

Lemma 3.1. For any environment u and any k,n>0, 

P^Xn = k)= ip^(Tfc+i = n + l). 

Proof. Note that Puj{Tk > m) = uj^^Pu}{Tk = m + 1) for all m > 0. Because 
Tk and r^, are independent, we find by applying (3.1) and conditioning on 
Tk that 

PuiXn = k) = Puj(Tk <n,Tk+Tk> n) 

= E^l{T^:<n}Puj{Tk +Tk> n\Tk) 

= UJ'^^E^l{T^:<n}Piu{Tk +Tk=n+l\Tk) 

= uj^^P^{Tk <n,Tk + Tk=n + 1). 

This implies that claim, because Tk+i = Tk+Tk and > 1 almost surely. □ 
3.2 Proof of Theorem 2.3 

By (2.4) we see that \ai^{Tk) < ck"^^ for some constants c > and A < 
1/2, so that X^fcLo ^^^'^('^'s) ^ "^tiich implies as a consequence of 
Kolmogorov's variance criterion (e.g. [23, Cor. 4.22]) that k~^ 'Y^^Z^i'^j ~ 
ojj^) — )• almost surely. By (2.5), we conclude that the hitting times of the 
walk satisfy the following law of large numbers: 

k~^Tk —7- /i almost surely. (3-3) 

Further, by (3.1), we see that Xn = T^{n + 1) — 1, where T"*" denotes the 
generalized inverse of the sequence {Tk) defined in Appendix A. Lemma A.l 
together with (3.3) now shows that n~'^T^{n) — t- /x^-*^ almost surely, and 
the proof of Theorem 2.3 is complete. □ 
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3.3 Outline of proof for the local limit theorem 



The starting point of the proof is Lemma 3.1 in Section 3.1, which reduces 
the problem into analyzing the probability mass function of the hitting times 
Tfc. Because is a sum of independent random variables, we may apply 
a classical local limit theorem (e.g. Petrov [""']) in Section 3.4 to conclude 
that PojiTk = n) ^ fk{n) for large values of k, where 



= (2w^)-V^ e , (3.4) 

and fik and denote the mean and the variance of given by (3.2). 

In the rest of the proof we need to transform the Gaussian density fk 
of the time variable into a Gaussian density of the space variable. This will 
be accomplished in two steps. We show in Section 3.5 that fk{n) gk{n), 
where 

gk{n) = (27rncjV/")"^/^ e ^"-'/'^ , (3.5) 
and further in Section 3.6 that gk{n) ~ fi^^hn{k), where 

(fc-fcn)^ 

hn{k) = (27rncrV/U^)"^/^ e ^^^^^ (3.6) 

is the Gaussian density appearing in the statement of Theorem 2.4 (we write 
kn in place of A;^ for convenience). 

The approximation gk{n) w fi^^hn{k) is a subtle part in the argument, 
where we have been guided by the following intuition. By approximating 
fikn ~ n, we obtain iXk - n ■^ii jik - IJ'k^ = Yleltl'^ where j = k - kn- 
Therefore, the difference of the square roots of the exponents in (3.5) and 
(3.6) can be approximated by 



(/ifc - n) iJ,{k - k, 



n ) 



(2naV^)-^/' (^7'-/^)- 

£=hn 



Condition (2.8) has been tailored to guarantee that the above difference is 
small for values of k such that iij. ^ n. 

An essential technical complication is the fact that the modulating factor 
uj~^^ in (2.10) (see also Lemma 3.1) is unbounded. To overcome this difficulty, 
we have strived to obtain sharp estimates. Identifying the set of reasonable 
sufficient assumptions (2.4) - (2.8) has played a crucial role in the proof. 



3.4 Local limit theorem for the hitting times 

By applying a classical local limit theorem for sums of independent random 
variables (Petrov [ Theorem VII. 5]; see also Davis and McDonald [13]), 
we obtain the following local limit theorem for the hitting times. 
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Lemma 3.2. For any environment uj satisfying (2.6) and (2.7), 

sup \P^{Tk = n)- fk{n)\ = 0{k-^), (3.7) 

n>0 

where the functions fk are defined by (3.4). 
Proof. Observe first that by (2.6), 

liminf A;"-^ Var(^(Tj) = ct^ > 0, 

j<k 

and that 

hmsup k"^ E^\Tj — E^Tjf < oo 

because of E^\Tj — E^Tj\^ < SooJ^ and (2.7). Therefore, to apply [-31, The- 
orem VII. 5], it suffices to verify (note that Pu}{Tj = 1) > Pu}{Tj = m) for all 
j and m) that 

1^ E^-(^^- = = 2) ^ oo. (3.8) 

By writing (^^^(1 — ^j) = xjyj, where xj = w^^^(l — ^j)^^^ yj = 
Uj ^"^^^ (l—ujj)^/^ , and applying Holder's inequality with conjugate exponents 
5 and 5/4, we see that 

Y,u^f{l-u^,)< Iy.^]{1-uA K;a;-3(l-a;,) 
j<k \j<k J \j<k 

After dividing both sides above by k, and applying (2.6) and (2.7), we see 
that 

liminf A;"-^ ^tj|(l - Uj) > 0, 

j<k 

which implies (3.8), because Pij{Tj = l)Puj{Tj = 2) = Ci;|(l — loj). The 

claim now follows by applying [ =1 , Theorem VII. 5] and recalling that o"^^ = 
0(fc-i/2) by (2.6). □ 

The following result is an analogue of Lemma 3.2, where the time variable 
n instead of the space variable k tends to infinity. 

Lemma 3.3. For any A G [0,1/2) and any environment uj satisfying (2.5), 
(2.6), and (2.7), the Gaussian densities fk defined by (3.4) satisfy 



snpk^\PUTk = n) - /fc(n)| = 0(n 

k>l 



A-l^ 
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Proof. Fix an e G (0,^ and note that 



sup e\PUTk = n)- fk{n)\ = 0^-^) 

k>en 



(3.9) 



by Lemma 3.2. Therefore, we only need to analyze Puj{Tk = n) and for 
1 < A; < en. As a prelimmary, note that, because e < fi~^ and — )• ^ by 
(2.5), we may fix a constant ci < 1 and an integer no such that fi[en] ^ cin 
for all n > uq. 

Assume now that k < en and n > ng. Then n — /i^ > n — > (1 — ci)n, 
and moreover, o"^ < u^en, where o"^ = sup£>;^((T^/^) is finite by (2.6). 
Therefore, 

{n- fikf/al>C2n, (3.10) 

where C2 = (1— ci)^/ (cr^e) > 0. A rough estimate together with Chebyshev's 
inequality shows that 

Puj{Tk =n) < P^{\Tk - lJLk\ >\n- //fc|) < (n - iik)~'^crl, 
so by (3.10), we conclude that 

k^P^{Tk = n) < c^^€\^~\ (3.11) 
By applying (3.10) once more, we see that 

k'^ fk{n) < (27ra2 )-i/2fcA-i/2g-c2n/2 < (2vra2 )-i/2e-c2n/2^ (3 -^2) 

where o"^ = inf£>i(cj|/£) is strictly positive by (2.6). The proof is now 
completed by combining the estimates (3.11) and (3.12) with (3.9). □ 

3.5 Variance of the hitting times 

Lemma 3.4. For any A G [0, 1/2) and any environment oj satisfying (2.5) 
and (2.6), 

max k^\fk{n) - gk{n)\ = o{n-^l^), (3.13) 

l<fc<n 

where the functions fk and are defined by (3.4) and (3.5), respectively. 

Proof. The proof is split into two parts according to whether or not |n— < 
ti(n log n)^/^, where li is a large constant to be determined later. 

(i) Assume that A; < n is such that |n — < u(nlogn)^/^. Using the 
triangle inequality and the inequality |e~^ — | < |x — y|, we see that 

(27rn)V2|/,(n)-5,(n)| 



< 



< 1 + {2a^/i^) 



-1/2 I"- ~ l^k\ 
n 
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Consequently, assuming that n is large enough so that 'u(logn)^/^ > 1 
(2^n)i/2 \fk{n)-gk{n)\ < ciu{logn)'/' \{n/cjl)'/' - {pi/a^f^ 

where ci = 1 + {2(j'^ / fiY'^/'^ . 

Observe next that, assuming n is large enough so that n — u(nlogn)^/^ > 
n/2, 

k = {fik/ky\n - (n - fXk)) > {2fi+y'n, (3.14) 
where /i+ = supgyi{iJ,i/i). Further, 

Now using (3.14) we find that al > c^'^n for C2 = 2/i+/(7^. Therefore, 

|n/a2 - /x/a^l = |(n - ^fe)/ai + d{k)\ < cstzn" (log n)!/^ + |d(fc)|, 
where 

As a consequence, 



n^/^A;^ |/fc(n) -5fc(n)| < C3un^(log n)^/^ (^nn-^/^^^Qg ^^1/2 ^ |^(^)| 

(3.15) 

for all large enough n, where C3 = cimax(c2,l). By combining (2.5) and 
(2.6), we find that \d{k)\ = o{k-^ (log k)''^/^). Because (2^+)-in < k <n 
by (3.14), we conclude that the right side above tends to zero as n — )• 00. 

(ii) Assume now that k < n is such that |n — > n(n log n)^/^. Note 
that cr'^k < al < (T\k, where o"^ = m.ii>i{a'j/t) and a\ = sup£>;^(cj|/^) are 
finite and strictly positive by (2.6). Therefore, the exponent in the definition 
of fk is bounded by 



la} 



> (2o-^)-^n2logn 



'k 

Consequently, 

I'^^k^ fkin) < 04^:^-^/2^1/2-^5^.^ < c4ni/2-'=5«^ , (3.I6) 



n 



where C4 = (27r cr^) ^/"^ and C5 = (2cj^) ^. For the function we immedi- 
ately see that 

n'/^k'^gkin) < c^n^l''-''^\ (3.17) 

where cg = (2'Ka'^ j and C7 = {2a'^ / p)~^ . The proof is now finished 
by choosing u > large enough so that c^v? > 1/2 and c-ju^ > 1/2, and 
combining the estimates (3.16) and (3.17) with (3.15). □ 
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3.6 From hitting times to walk locations 



Lemma 3.5. For any environment uo satisfying (2.5) and (2.8) for some 
AG [0,1/2), 

max - /i ^hn{k)\=o{n ^/^), 

l<fc<n 

where the functions and /i„ are defined by (3.4) and (3.6), respectively. 

Proof. We will show the claim by treating separately the cases I | ^ 
ub{kn) and \k — kn\ > ub{kn), where b{k) = (klogk)^/"^ is as in (2.8), and 
It > is a large constant to be determined later. 

(i) Assume that k < nis such that \k—kn\ < ub{kn). Using the inequality 
\e~^ — e~y I < |x — y|, we see that 

|gfc(n) - /i~^/i„(/c)| < cin~'^\d{k,n)\, 

where ci = 7r~^/'^{2a'^ / n)^^ , and d{k,n) = (n — fij.) — ^{kn — k). Note that 



\d{k, n)\ 



hfi 1 



l=k 



< 



kji 1 



£=k 



where the latter inequality is due to fik^ — ^i^^^i < n < ^1^^. To analyze the 
right side of (3.18), observe that, because kn — )• oo, we see using (2.8) that 



max 

k:\k— kn\<ub{kn) 



kn 1 



i=k 



max 

j:\j\<ub{k„) 



kn+j-l 

E 

i. — hfi 



-1 



1/2-An 



The above limiting relation also shows (substitute j = 1) that = 
o{k}/'^ Because kn/n — )• fi~^^ (by (2.5) and Lemma A.l), it follows that 

\d{k,n)\=o{n'/^-^), 



max 

k:\k—kn\<ub{kn) 



and therefore, 



,1/2 



n"'" max k^\gk{n) — fj, ^hn{k)\ 

l<k<n:\k—kn\<ub{kn) 



0. 



(3.19) 



(ii) Assume that A; < n is such that \k — kn\ > ub{kn) for some u > 2 
and n > no, where uq has been chosen large enough so that kn > \^~^n 
for all n > no (this is possible by virtue of (2.5) and Lemma A.l). Observe 
first that, because uj^^ > 1 and \k — kn\ > 2, we see by Lemma A. 2 that 
\l^k — n] > |A; — A:„| — 1 > ^\k — kn\. Hence, the exponent in the definition of 
gkin) is bounded by 

A^'t ) ^Ti) 2ij 2i/N 2i 

> Cl ; > ClU log kn > —ClU log(2/i) + ClU iogU, 



kry 
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where ci = (16cj^) ^. The same bound is also vahd for the exponent in the 
definition of hn{k), because fi^ > ^. Therefore, we obtain 

n'/^k^\gk{n) - ii-^K{k)\ < C2(2/x)^i"'n^-'=i"', (3.20) 

where C2 = 2(27r(T^//i)~^/^. The right side above tends to zero as n — )• oo, if 
in addition to u > 2, we also require that u > (A/ci)^/^. The proof is now 
completed by combining (3.19) and (3.20). □ 

Lemma 3.6. For any A > and any environment uo satisfying (2.5), there 
exists a constant c > such that the functions hn{k) defined by (3.6) satisfy 

snpk^hnik) = o(e-""). 

k>n 

Proof. Recalling (by (2.5) and Lemma A.l) that kn/n ^ < 1, we may 
fix a positive constant ci < 1 and an integer hq such that kn < cin for all 
n > uq. Assume now that k > n and n > uq. Then, k — kn > {1 — ci)k, so 
that the exponent in the expression of hn{k) is bounded by 

(fe_fc„)2 (1-ci)^ 
2na7/i3 - 2a2^3 

As a consequence, 

where C2 = {27ra'^ / fx^y^^'^ and C3 = (1 — ci)'^ {2a'^ / fi^)^^ . Because the 
function t 1— )• t^e~'^^^ is decreasing on the interval [A/c3,oo), we conclude 
that 

SW^k^hn{k) < C2n^"l/2g-C3n^ 
k>n 

for all n > max(no, A/cs), so the claim follows. □ 

3.7 Proof of Theorem 2.4 

By combining Lemmas 3.3, 3.4, and 3.5 we see that 

max k^ |-Pcj(Tfc = n) — fi~^hn{k)\ = o(n~"^/^), 

l<fc<n 

where the functions /i„ are defined by (3.6). Hence by Lemma 3.6 and the 
fact that Puj{Tk = n) = for all k > n, we conclude that 

supA:^ \PUTk = n) - fi-^hn{k)\ = o(n-V2). 

k>l 

Further, because < cik"^ < ci{k + 1)'^ by (2.4), we see by Lemma 3.1 
that 

sup \PUXn = k)- {uj-y ^i)hn+i{k + 1))| = o(n-i/2). 

fc>0 
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We complete the proof of Theorem 2.4 by showing below that 

sup/i:^|/i„+i(A; + 1) - /i„(A:)| = o(n"^/2)_ (3 21) 



fc>0 



Let us write hn{k) = cin ^^^e '^^"fc.", where ci = {2-iTa'^ / fi^) C2 



(2(T7/i^) \ and a^. 



n 



'^/'^{k - kn). Note that 



ak+i,n+i-ak,n = {n + l)-^l'^{kn - kn+i + l) + {k- kn){{n + l)-^/^ -n~^/^). 

Note that < kn+i — fc„ < 1 (Lemma A. 2) and \k — kn\ < (1 + C3)n for 
all k < n, where C3 = sup^>o(A:£/£) is finite by Lemma A.l. Therefore by 
applying the inequality n~^/^ — (n + < ^n~^^'^, we see that 

|afc+i,„+i - ak,n\ < (2 + (1 + C3)/2)n-^/2 



for all < n. This estimate combined with the inequality \e 
\x — y\ now shows that 



maxk^\hn+i{k + 1) - /i„(A;)| = 0(n^~^). 

k<n 

Together with Lemma 3.6, we now conclude the validity of (3.21), and the 
proof of Theorem 2.4 is complete. □ 

3.8 Proof of Theorem 2.5 

It suffices to show that 



for all X < y, where o"^ = j [j? . Note that the left side above can be written 
as X^jtg/^ -Pa;(^n = fc); where the set /„ = Inix,y) is defined by 

In = {k ^ : \xa\fn\ + 1 < k — kn < [u^V^} > 

and by Theorem 2.4, 

PUXn = k)- ^k^^nik) ^ 0. 

Note that for any real numbers x < y and any sequence {ak)kez, 



a[t] dt- Y 

k=\x]+l 



< lar^il + la 



[y]\ 



(3.22) 
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By applying (3.22), using (2.4), and performing a change of variables, we 
see that 

u:~^K{k) = r ^' u-^^^^^hniK + M) du + 0(n^-i/2) 



k&In 



y 

f^Vn I i^ll^^^^^K{kn + [cfty/n]) dt + o(l). 



where 0(n^"^/2) = o(l) due to A < 1/2. Further, - [tf\ < 2\t\ shows 
that 

KiK + [atV^] = (2w2n)^V2exp (-^^^ 
= (27ra2n)-V2e-*V2o(eM/(Sv^)) 

because we assume that t G ix,y). 

The next lemma finishes the proof by showing that 

' ''\-^'^..^/^^)e-''^'dt^^^e-^y^dt. 



Jx ^•"+[^*v^'^^ ^ 

The intuition is that, restricting integration to sufficiently small subinter- 
vals of (x, y) and making n large, e~* is virtually constant, while the 
fluctuations of uj~ . ,-, ^ average out in the integral. 

Lemma 3.7. On the interval {x,y), the probability measures defined by 
Tan{dt) = Z^"^ n^^uj^ dt converge weakly to the uniform probability 

measure m(dt) = {y — x)~^ dt and the normalizing factor Zn ^ y — x. 

Proof. We must show that m„((x,s]) — t- m((x,s]) for all s € {x,y). By a 
change of variables and (3.22) we see that 

because uj'^^ = 0{k^) and A < ^. The center of In{x,s) is c„ = /i^^n(l + 
o(l)). For u > large enough and b{n) as in (2.8), |/„(x,s)| = [sa^/n] — 
[x^^^/n\ < ub{n). Therefore, (2.8) implies that 



n 



^1/2 



^k^=n ^/^\In{x,s)\fi + oil). 



Because n ^/^|I,i(x, s)| — )■ (s — x)(5", we get /i ^ uj^^_^^^^^-^dt ^ s — x, and 
Zn ^ y — X follows by taking s = y. This finishes the proof. □ 
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4 Proofs for quenched random environments 



In this section we prove Theorem 2.7 for the random walk (Xn) in a quenched 
random environment satisfying the assumptions (2.13)-(2.14). Section 4.1 
contains some preUminary facts on the growth rate and moving averages of 
stationary sequences, and the proof of Theorem 2.7 is given in Section 4.2. 



4.1 Growth rate and moving averages of stationary sequences 

The following result establishes a bound on the growth rate of a stationary 
sequence in terms of its moments. 

Lemma 4.1. Let {(,o,(,i, ■ ■ ■) be a stationary random sequence such that 
El^ol'' < CO for some q > 0. Then = o{k^^'^) almost surely. 

Proof. Fix e > 0. By stationarity and Fubini's theorem we see that 

oo oo 

Y,nk-'^'\Ck\ > e) = ^P((e-i|Co|)^ > A:) = E[{e-'\Co\y] < oo, 

fc=0 fe=0 

where [x] denotes the integer part of x. Because e was arbitrarily chosen, 
the claim follows by the Borel-Cantelli lemma. □ 

The next result analyzes the moving averages of a stationary random 
sequence in terms of its moments and mixing rate. The result, proven with 
the help Peligrad's law of large numbers [29] (see also Bingham [ J] for a 
nice survey), is similar in spirit to Kiesel [26, Thm 1], but tailored to fit our 
needs. 

Lemma 4.2. Let {(,o,(,i, ■ ■ ■) be a stationary random sequence such that 
El^ol"^ < CO for some q > I, and for which the mixing coefficients as defined 
in (2.15) satisfy 

oo 

^ 0^/''(2") < oo for some k>2. (4.1) 
n=l 

Then 



max 



k+j 
'=k+l 



(4.2) 



almost surely for any s > 0, u > 0, and 



, 1 + s s 1 , , , , 

.>max(--, ^ + ^). (4.3) 
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Proof. Denote the left side of (4.2) by M^, and fix an arbitrary e > 0. Note 
that by stationarity, 



^P(fc-M, > e) = ElP(i<-<-,. I^.l > (4-4) 

k=l k=l 

where Sj = J2e=ii^i ~ ^Co)- The right side in (4.4) is finite if and only if 



max IS",! > et'^ ] dt < oo, 



and by a change of variables, the above is further equivalent to 
^ ti/«-ip(^ max^ \Sj\ > eif dt < oo. 



where ei = u ^/*e. We conclude that Yl'kLi ^{^^ > e) is finite if and 

only if 

oo 

J]]n"P-2p(^ max \Sj\ > ein") < oo, (4.5) 

n=l — J— 

where a = r/s and p = (1 + s)/r. 

The parameters r and s have been transformed into a and p to conform 
with the notations used in ]. Observe that ap > 1, and moreover, (4.3) 
implies that p < q, a > 1/2, and [(ap — l)/(a — 1/2)] + 1 < k. Hence, 
Peligrad's law of large numbers Thm 2] yields the validity of (4.5), and 
therefore, Yl'k^i ^{k~'^Mk > e) is finite. Because e was arbitrarily chosen, 
the claim follows by the Borel-Cantelli lemma. □ 



4.2 Proof of Theorem 2.7 

By virtue of Theorems 2.3-2.5, we only need to verify that the regularity 
conditions (2.4) - (2.8) are valid for P-almost every realization of the random 
environment. 

Recall that Kujq < oo for some q > 5hy (2.13). Consequently, Lemma 4.1 
implies that = o{k^/'^), so that especially, (2.4) holds with X = 1/q. Be- 
cause q > 5, the random variables = oj^^ and s| = (1 — uJk)(jJ^'^ have 
finite second moments, so the assumption (2.14) on the mixing coefficients 
implies that the stationary random sequences (m^) and (s|) satisfy the law 
of the iterated logarithm [5, Section 12]. As a consequence, (2.5) and (2.6) 
are valid P-almost surely with n = Emg > 1 and = Esq > 0. 

Because Ewq^ is finite, the pointwise ergodic theorem implies that 
=o^k^ converges, so that (2.7) holds. 

To verify (2.8), note that because (p is decreasing, the sum in (4.1) equals 
joo^i/^^2N)dx < /;^<^i/'^([2^-i])dx = {log2)-^ (l)^/^{[y])y-Uy. Us- 
ing Holder's inequality, the last integral can be bounded from above by 
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(/r (t^'^'^M) dyf" (/r y'^dyf" withp = (1-2/k)-i, which is finite due 



to (2.14). Therefore, condition (4.1) of Lemma 4.2 holds for any k> 2. Next, 
because g > 5, we may choose an exponent s > 1/2 such that s < 1 — 2q~^ 
and s < q/2 — 2. Our choice of s imphes that the exponent r = 1/2 — 
satisfies (4.3) for some large enough k. Hence by Lemma 4.2, it follows that 



5 Conclusions 

We studied the propagation of a particle in a one-dimensional inhomoge- 
neous medium, where the motion is induced by chaotic and fully determin- 
istic local rules, and the initial condition is the sole source of randomness. 
The spatially varying local rules constitute an environment, which is frozen 
during the particle's lifetime. This model falls into the framework of ex- 
tended dynamical systems which lack physically observable invariant mea- 
sures. Defining the local rules via piecewise affine maps allows to reduce the 
model to a simple unidirectional random walk on the integers. 

The main result of the paper shows that the probability mass function 
of the random walk approaches a modulated Gaussian density, where the 
modulating factor is explicitly given in terms of the local properties of the 
environment. In contrast, when looking at the walk over a coarser diffu- 
sive space scale, the non-Gaussian modulating factor averages out asymp- 
totically, and the distribution of the walk approaches a standard Gaussian 
distribution. 

Although our analysis is restricted to a special instance of a random walk, 
we believe that the obtained results could serve as useful benchmarks when 
testing hypotheses concerning more general extended dynamical systems and 
random walks in random environments. 

A A generalized inverse 

Let (a(/c))^Q be an increasing sequence such that a(0) = and 
limfc_^oo a(fc) = oo, and define its generalized inverse by 



Then also a^(0) = and lim„_>oo o"^('i-) = oo. The following result sum- 
marizes some basic properties of the inverse. 




almost surely for all u > 0. Because (fclog k)^!'^ 
holds with A = 1/g. 



o(A;^), it follows that (2.8) 

□ 




(A.l) 
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Lemma A.l. If \\mk^o^ a{k) / k = fj,, then lim„_j.oo o^(f^)/'^ = 1/^- 

Proof. By definition, a{a^{n) — 1) < n < a(a^(n)) for all n such that 
a*^{n) > 0, so that 

a(a^(n) — 1) n ^ a{a^{n)) 

Because a^(n) — ?■ oo, it follows that a{a^ (n)) / (n) — )• /x. Therefore, the 
above bounds imply n/a^{n) — )• /i, and the proof is complete. □ 

Lemma A. 2. Assume that c = infjt(a(fc + 1) — a{k)) > 0. Then for all 
positive integers k and n, 

\a{k) -n\>c{\k- a^{n)\-l), (A.2) 

and 

a^(n + 1) - a^(n) < [c-^]. (A.3) 

Proof. If A; > a^(n), then a{k) — n> a{k) — a{a^{n)) > c{k — a'^{n)). If 
k < a'*~(n), then n — a{k) > a(a^(n) — 1) — a{k) > c{a^[n) — k — 1). Hence 
(A.2) follows. 

To prove (A.3), fix a positive integer n, and denote H. = a^[n). Then 
a{t) > n and a{£ + [c"^] ) - a{e) > c\c-^] > 1, so that a{£ + \c~^] ) > n + 1. 
Hence, a^(n + 1) <£+ [c"^]. □ 
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